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ABSTRACT 

We add the a— effect in the dynamo driven accretion disk model 
proposed by Tout & Pringle (1992), i.e., a dynamo model depends on the 
physical processes such as Parker instability, Balbus-Hawley instability, 
magnetic field reconnection and a — uj mean field dynamo as well. The 
a— effect in the dynamo mechanism is determined by the strength of 
turbulence of the accretion flow. When the turbulent Mach number Mt is 
less than 0.25, the solutions of the magnetic fields oscillate around their 
equilibrium values. The increase of the value of Mt makes the amplitude 
of the oscillation smaller and the period longer, but does not affect the 
equilibrium values. The Shakura-Sunyaev viscosity parameter ags oscillates 
around the equilibrium value of 0.33. When the turbulent Mach number 
Mt is larger than 0.25, the magnetic field components reach a stable state. 
In the non-linear dynamo region, the critical turbulent Mach number Mt 
is 0.44 rather than 0.25. The oscillating magnetic fields and viscosity 
parameter can explain the basic properties of the dwarf nova eruptions and 
some properties of quiescent disks (Armitage et al. 1996). 
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turbulence 
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1. 



Introduction 



There is observational evidence for the presence of magnetic fields in accretion 
disks, and they influence the dynamics of the system, such as angular momentum 
transport and poloidal outflow of the disk gas. The nature and magnitude of the 
viscosity is uncertain in the current theoretical descriptions of standard accretion disk, 
but almost all detailed modeling on the structure and evolution of accretion disks 
depends on the value of the viscosity (Duschl et al. 1997). In general, the molecular 
viscosity is inadequate and some kind of turbulent viscosity is required. Shakura & 
Sunyaev (1973) introduced a single dimensionless parameter which we denote here 
as ass to describe all the unknown physics about the viscosity, so that the kinetic 
viscosity u is written as 



where is the sound speed, Q the angular velocity in the disk. 
When a magnetized accretion disk is concerned, the dimensionless measure of the 
strength of the viscosity ass, i-e., the Shakura-Sunyaev viscosity, is determined by the 
magnetic field in the form 



where Br and are the radial and azimuthal components of the magnetic fields, p the 
density of the disk, and Vai notes the Alfven speed corresponding i component of the 
magnetic field. Most investigators think that the magnetic fields in the accretion disk 
are substained by a magnetic dynamo. Tout & Pringle (1992) set forward a physical 
model simply considering the Balbus-Hawley (B-H) instabihty, Parker instability and 
Shearing rotation without a— effect. They claimed that only the three well-known 
physical processes are included in their model. This model was used to explain the basic 
properties of dwarf nova eruptions and some properties of quiescent disks (Armitage et 



u = assC'jn, 




ass = BrB^/A7:pCl = VArVA^/Cl 



(2) 
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al. 1996). But the a effect which arises from the correlation of small scale turbulenct 
velocity and magnetic fields is important in maintaining the dynamo action by relating 
the mean electrical current arising in helical turbulence to the mean magnetic field. It 
plays the key role in the current studies of astrophysical dynamo and many authors 
have studied it in detail (Krause & Radler 1980, Stepinski & Levy 1994, 1991, Elstner 
et al. 1994). Furthermore, the nonhnear effect of the a— quenchings were studied in 
detail (Schultz et al. 1994, Riidiger & Schultz 1997, Covas et al. 1997). ol-oj dynamo 
models have been shown to be capable of producing a number of important features 
of different astrophysical objects. The presence of strong differential rotation plus a 
vertical density gradient will give rise to the turbulence helicity. Considering that a— 
effect arising from the helicity of turbulence is commonly accepted in Keplerian disks 
as well as stellars (Krause & Radler 1980, Pudritz 1981a b, Mangalam & Subramanian 
1994, Vishniac & Brandenburgl997, Reyes- Ruiz & Stepinski 1997), we add a— effect 
in their model and reconsider the behaviours of magnetic fields and the viscosity 
parameter in this paper. Instead of Balbus-Hawlcy instability considered in this model, 
the magnetic buoyant combined with Coriolis twist was adopted to close the dynamo 
cycle by Rozyczka et al. (1995). They found that the magnetic fields reach a stable 
saturation state. We feel that their scenario is approciate for a magnetized disk in 
which the Balbus-Hawley instability is restained. 

In Section 2, we write down the full dynamo equations. The physical processes cover 
the Parker instability, Balbus-Hawley instability, a — io dynamo mechanism, and 
magnetic reconnection. The dynamo equations are solved in Section 3. Finally, we 
present conclusion and discussion in Section 4. 



2. The model 
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2.1. Parker instability 

Parker instability is a kind of Rayleigh- Taylor instability, i.e., interchange 
instability, introduced by magnetic buoyancy (Parker 1979). It was found that the 
gas disk supported in part by the magnetic field against vertical external gravity is 
unstable against the long wavelength disturbance along the field lines. The horizontal 
components of the magnetic field are denoted by Br and 5</,, and the vertical component 
by Bz- The Parker instability leads to a loss of horizontal magnetic field and converts 
it to a vertical field. If the horizontal magnetic field is dominated by the azimuthal 
component 5^, we have 

dBz B(h dBj. Br dBj, Ba, 



dt Tp ' dt Tp^ dt Tp 
The growth time of the Parker instabihty is 



(3) 



Tp = r]H/VA4>, (4) 

where H is the half-thickness of the disk and r; ~ 3 in a non-selfgravitating accretion 
disk according to the calculation given by Houriuchi et al. (1988). The wavelength of 
the instability in the azimuthal direction is 

Ap0 = e^, (5) 

where ^ ~ 8. 



2.2. Balbus-Hawley instability 

In recent years, Balbus & Hawley et al. published a series papers concerning the 
stabihty of weakly magnetized disks (Balbus & Hawley 1991a b, 1992). It was found 

that the shear instability is local and extremely powerful provided that the field energy 
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density is much less than the thermal energy density. The maximal growth rate is given 
by the local Oort A-value of the disk, which is not only independent of the magnetic 
field strength but independent of field geometry as well (Balbus & Hawley 1992). If 
a fluid element is outwardly displaced in a differentially rotating disk threaded by a 
vertical magnetic field, rigid rotation will be enforced because it is elastically tethered 
by a magnetic field, i.e., the field is trying to force the element to rotate too fast for 
its new radial location if ^ < 0. The excess centrifugal force drives the element still 
farther outward. Therefore, the effect of the instability is to tap the energy present in 
the shear ffow and to use it to generate radial field from the initial vertical field, i.e., 
the vertical magnetic field results in the amplification of the radial component of 
the magnetic field Bj. because of Balbus-Hawley instability 

dBr 



dt 



(6) 



'Jmax 

10, 



(73-1)2 



1/2 



where is the growth rate of the instability, 

Vaz/Cs < \/2/7r, 
, ^/tt < < \/6/7r, 

VaJCs > v^/tt. 



(7) 



here 7„ 



\Q, which equals the Oort A-value of a Keplerian disk. The scale of the 



instability in the vertical direction, Xbh, is given by 



Xbh 




Va. 

Cs 


> 


V2 


2n 


[ V2Cs 


Cs 


< 


TT 



(8) 



We can see that Balbus-Hawley instabihty is cut off if the magnetic field is strong 
enough. 
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2.3. The a — uj dynamo 



The mean field dynamo theory has been successful in many kinds of astrophysical 
objects, and the a effect which arises from the correlation of small scale turbulenct 
velocity and magnetic fields is important in maintaining the dynamo action by relating 
the mean electrical current arising in helical turbulence to the mean magnetic field. It 
plays the key role in the current studies of astrophysical dynamo (Ma & Wang 1995, 
Stepinski & Levy 1991, Krause & Radler 1980). The a — u dynamo models have 
been shown to be capable of producing a number of important features of different 
astrophysical objects. The presence of strong differential rotation plus a vertical density 
gradient will give rise to the turbulence helicity. The radial component of the magnetic 
field is amplified by means of the a— effect following the equation 



The parameter a denotes the product of the mean helicity and the correlation time of 
the turbulent flows. In term of the disk's angular velocity and turbulent Mach number, 
we take the form (Stepinski & Levy 1991) 



where Mt = Vf/Cg is the Mach number of the turbulence and Vt the turbulent velocity. 
The recent numerical simulation indicated that the nonmagnetized astrophyical 
accretion disks are both linear and nonlinearly stable to shearing instabilities, thus 
ruled out any kind of self-generated hydrodynamical turbulence (Balbus et al. 1996). If 
the source of the turbulence is the Balbus- Hawley instability in a weakly magnetized 
disk, the turbulent Mach numbers fall in the region between 0.1 and 0.25 (Hawley 
et al. 1996). So, the turbulent Mach number Mt is adopted the order of 0.1 in our 
calculations. 



dBr 




(9) 



dt 



a = HM^n, 



(10) 
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At the same time, the amphfication of azimuthal magnetic field results from the 
a;— effect, i.e., the differential rotation, satisfying the equation 

dBs dn^ 3^^ 

-^^r^B, = -aB., (11) 
Finally, we have closed the cycle for the dynamo to work. 



2.4. Dissipation of the magnetic energy 

The dominant flux loss mechanism is the reconnection of the vertical component 
Bz in the radial direction. Consider two patches of of opposite sign coming together 
and reconnecting within distance of A^ec- We have a term of dissipation in the vertical 
equation in the following form (Tout & Pringle 1992): 

The magnetic flux is removed from the disk owing to reconnection at a rate 

- = - + ^. (13) 

' rec ' rec ' rec 

where r^gc is determined by the length scales of Parker instabihty, B-H instabihty and 
the shearing of the rotation, 

while r^gg results from the shearing itself (Tout & Pringle 1992), 

- {^) ■ <^^' 

Here F"^ ~ ln(3fJm), where 3?^ is the magnetic Reynolds number and F is expected to 
be in the range 0.01 to 0.1 (Tout & Pringle 1992). According to the analysis made by 
Tout & Pringle, we adopt Xrec(i> — 0-5Xp(I>Xbh/H. 
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2.5. The full dynamo equations 



Concluding all the physical processes stated above, and defining a dimensionless 
time T = -\/2^/^,7 — 7/^ and a dimensionless velocity Wi — Vai/Cs 
iVAi — ^il s/^T^Pi i — T^,<l>,z), we get the full dynamo equations in the spatially 
local approximation 

1/2 



dWr 



. 1 _ A-w^ v 



(16) 



dws 3\/2 2 



(17) 



> V2/ 



TT. 



(18) 



We note that all the spatial effects are ignored in these equations. As compared with 
Tout & Pringle's results, the terms including Mf present the a effect in Eq.(16), which 
are not invoked in their model. In the next section, we will analyse the solutions of the 
full dynamo equations. 

3. Solutions of the dynamo equations 



3.1. The equilibrium solutions 



When the turbulence is weak, the process is dominated by B-H instability and 
Parker instability, so the equilibrium analysis should be the same as Tout & Pringle's 
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results. On the other hand, if the turbulence is so strong that Mi > 0.25, we will 

see from the numerical calculation that the a— effect takes over, and the equilibrium 
solutions read as 

^ = V2vMl (19) 

^ = VSrjM,, (20) 

1^ = 3l/42-3/4^1/2p-l/2^1/2^3/2_ ^21) 

In this regime, we estimate the viscosity according Eq.(2) 

ass = Vev^Mf. (22) 



3.2. The numerical solution 

Based on the statements in Sec. 2, we take the parameters F = 0.1,77 = 3,^ = 
8,lmax — 0.75 and the initial condition i/;,. = = = 0.01 at r in our 
calculations. When the turbulent Mach numbers are set to be 0.1, 0.2, 0.24 and 0.30, 
the time dependence of the magnetic field components are shown in Fig.l - Fig.4 
respectively. We can see from the figures that the dynamos are dominated by B-H 
instability, Parker instability and the shearing motion of the disk if the turbulent 
Mach Numbers are less than 0.25, and the magnetic fields oscillate around their 
equilibrium values. The Mach number of the turbulence does not affect the equilibrium 
value of the amplified magnetic fields so much, but the amplitudes of the oscillation 
of the amplified magnetic field decrease meanwhile the periods increase as the Mach 
number increases. If the turbulent Mach number is larger than 0.25, the a — ou 
dynamo takes over and the oscillation of the amplified magnetic field disappears, 
and the equihbrium values depend on Mt as Eq. (21) — (23). The corresponding 
time-dependent behaviors of the viscosity parameter ass are shown in Fig. 5, where the 
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turbulent Mach numbers are set to be 0.1, 0.2 and 0.25 respectively. Wc can sec that 
the time averaged value of ass is about 0.33, which does not depend on the turbulent 
Mach number Mf. 

3.3. The non-linecir dynamo 

When the back-action of the amplified magnetic field on the turbulent motion of 
the fluids is considered, the mean helicity of the turbulent flow becomes 

a' = q;*(S2), (23) 

where the function "^{B"^) = (1 + V^/C^)"^ is quadratic in magnetic field, representing 
the so called a— quenching(Schultz et al. 1994, Riidigcr ct al. 1997, Covas et al. 1997), 
and a is determined by Eq. (10). The critical turbulent Mach number increases to 
Mt ~ 0.44, as compared with the hnear model, it is 0.25. 

4. Conclusion and discussion 

Tout & Pringle (1992) put forward a physical model including three processes: 
the Parker instability, the Balbus-Hawley instability and magnetic field rcconnection, 
to study the magnetic field configurations and viscosity parameter. The magnetic 
dynamo origin for the viscosity based on this model was used to explain the eruptions 
of accretion disks and dwarf nova (Armitagc ct al. 1996). Considering the a effect has 
been well studied too, we adding the a effect in their dynamo model and obtain the 
results as follows. 

If the turbulent Mach number < 0.25 in the linear dynamo model, dynamos are 
dominated by the Balbus-Hawley instability and the Parker instability. When the 
dynamo works in this region, we have the time- averaged viscosity ass ~ 0.33. The 
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turbulent Mach number associated with the a — uj dynamo only affects the amplitude 
and the period of the oscillation of the amplified magnetic fields. The kinematic 
viscosity obtained from this simple model can be used to explain the outburst and 
quiescent phases of the disk around the Dwarf Nova (Tout & Pringle 1995, Armitage 
et al. 1996). The disk remains in a quiescent state with strong fields and low viscosity 
until sufficient mass has been added to restart the B-H instability. On the other 
hand, Gammie and Menou (1997) proposed an other scenario for the origin of episdic 
accretion in darf novae, in which the outburst cycle purely results from the grobal 
hydrodynamic instability and depends on the magnetic Reynolds number. This model 
is totally different from the standard disk instability model. 

But if the turbulent Mach number Mt > 0.25, the Balbus-Hawley instability is 
restrained and the oscillations of the magnetic field components disappear. That means 
that the equilibrium magnetic fields are determined by Mt alone. The corresponding 
viscosity is ass — V^rfMf. In this case, the disks should keep in a quiescent phase all 
the time, i.e., there would be no eruptions any more when the turbulence is so strong 
that the Mach number is larger than 0.25. But the turbulent Mach number depends 
on the temperature of the accretion flow. That means that a combination of the disk 
model and dynamo model is necessary to describe the behavior of the accretion disk 
more clearly. We note that the standard mean field dynamo theory does not cover the 
Balbus-Hawley instability. 

Therefore, if the turbulence in the disk is weak enough, the a— effect does not change 
the basic instability of the magnetized accretion disk described by Tout et al. (1992). 
The physical nature of turbulence in astrophysical objects is not clear, although most 
people believe that the turbulent Mach number in accretion disks or galaxies is of the 
order of 0.1 (Ruzmaikin et al. 1980, Stepinski et al. 1991, Moss et al. 1996). The three 
dimensional simulations of an aacretion disk indicate that the values of the turbulent 
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Mach number fall in the region between 0.1 and 0.25 (Table 4 of Hawley et al. 1996). 
Brandenburg et al. (1995) suggested that supersonic flows, initially generated by the 
Balbus-Hawley magnetic shear instability, regenerate a turbulent magnetic fleld, which, 
in turn, can reforce the turbulence. The flnal turbulence is not so strong based on 
their simulations either. So it seems that the turbulent dynamo is not so effective to 
suppress the Balbus-Hawley instability completely. In fact, when a nonlinear dynamo 
is concerned, the critical turbulent Mach number increases to be 0.44 rather than 0.25 
as compared to the linear model. We conclude that the a— effect does not change the 
basic conflguration of magnetic fields described by Tout & Pringle (1992). 
There are two facts in the model may need modification for further studies. The first 
is that the dissipation processes of magnetic fields in accretion disks is not clear to 
us, which are simplified to some degree in this model. The second is that the spatial 
behavior of a full magnetic dynamo model should be taken into account. 
We note that the vertical component of the magnetic field is important for the outfiow 
of jets in AGN (Yoshizawa & Yokoi 1993). Tout & Pringle (1996) demonstrated that 
a dynamo-generated field coupled with an inverse cascade process is able to produce 
sufficient field strength on large enough scales to drive a large-scale outfiow. The 
similar inverse cascade process of the magnetic fields from small scale to large scale was 
recently proposed by Kulsrud et al. (1997) to explain the origin for cosmic magnetic 
fields. 

The rotating magnetic fields in the accretion disk would introduce an electric field 
which can accelerate the charged particles. Some of the electrons will become run-away 
particles in the electronic field. The high energy electrons may emit X-rays or 7-rays 
from the disk. (Pustfl'nik & Ikhsanov 1994). 

GYM is supported by a fellowship of the Max-Planck Society. 
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Fig. 1. — Time dependence of the components of the magnetic fields through the dynamo 
cycle when Mt — 0.10. The solid line corresponds to the z-component Wz, the dashed 
line to Wtf, and the dotted line to Wr- 
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Fig. 2.— The same as Figure 1. but Mt = 0.20. 




Fig. 3.- 



The same as Figure 1. but Mj = 0.24. 
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Fig. 5. — The time-dependent behavior of the viscosity parameter ass- The dotted hne 
corresponds to Mt — 0.10, the dashed hne to Mt — 0.20 and the sohd hne to Mt — 0.25 



